Orbital ordering in undoped manganites using Jahn- Teller interaction 
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Abstract : To understand the orbital ordering of LaMnOz in the ground state, we study it in the A-type spin antiferro- 
magnetic state. We calculate the two-dimensional response functions associated with the Jahn- Teller Qi and Qz distortions 
and find that they diverge at the wavevector (7r,7r). Furthermore, the Q2 response function diverges faster which indicates that 
perhaps the ordering of the orbitals at low temperatures is dominated by the Q2 distortion. We calculate the ground state 
energy for the cases when only one of the Q2 and Q3 modes is excited cooperatively and find that the Q2 excited state yields 
lower energy. 
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1. Introduction 

Undoped manganites like LaMnO^ are the parent sys- 
tems for the colossal magnetoresistive materials. It is 
well known that orbital ordering occurs around 780 K 
resulting in a C-type orbital structure with two kinds 
of orbitals alternating on adjacent sites in the xy plane 
while like orbitals are stacked in the z direction [1]. As 
the temperature is further lowered to 140 K an A-type 
spin antiferromagnetic order sets in wherein the spins 
are ferromagnetically aligned in the xy plane with the 
spin coupling in the z direction being antiferromagnetic 
[2]. To explain the observed order several studies have 
been reported. These studies fall into two classes, with 
one class based on Jahn- Teller (JT) interaction [3-6] be- 
ing the main cause while the other class corresponds to 
Coulombic interaction [7-10] being the dominant cause 
for the observed order. However, it must be pointed out 
that the electron-phonon based studies do not point to a 
transparent mechanism. 

The purpose of the present paper is to use the frame- 
work developed earlier by one of us (Ref. [11]) and ana- 
lyze the orbital ordering by using a generalized Peierls 
instability approach. However, as compared to the 
one-dimensional Peierls charge density wave (CDW) ap- 
proach, our higher dimensional orbital density wave 
(ODW) analysis is more complicated on account of there 
being two e g orbitals (with inter-orbital hopping) and 
two response functions corresponding to the JT Q2 and 
Q3 distortions. We begin by assuming that the A-type 
antiferromagnetic ordering has already set in. Thus, 



on account of strong Hund's coupling, the transport 
is restricted to spin polarized electrons in two dimen- 
sions only. We assume that the Mn sites are fixed and 
that only the oxygen atoms between the Mn sites can 
move along the axis joining the neighboring Mn sites. 
We find that the non-interacting response functions X2,3 
corresponding to the Q2,3 modes diverge at wavevector 
Q = (jr, it). Thus one can expect local minima in en- 
ergy for states corresponding to ODW being excited at 
wavevector Q by the cooperative Q2 and Q3 modes. We 
find that ^(tTjI") diverges faster than Xsi 1 ^^) an d as 
expected the energy for the Q2 related ODW state has 
lower energy than the Q3 related ODW state. 

2. Mean-field approach to ODW 

We will now consider manganite systems with two e g 
orbitals per site and ignore spin. The Hamiltonian con- 
sists of the kinetic term, the ionic term, and the electron- 
ion interaction term. The kinetic term in momentum 
space is given by 
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where Bt = (b\ b\^) with b\ and 62 corresponding 
to the destruction operators for electrons with the or- 
thonormal wavefunctions ip x 2_ y 2 and V'3z 2 -r 2 respec- 
tively. Furthermore, T is a hermitian matrix with T^i = 
— 1.5t[cosp x +cospy], T 2 ,2 = —0M[cosp x +cosp y ], and 
TY2 = 0.5\/3t[cosp x — cos p y ]. The eigen values of the 



1 



kinetic energy are given by A^ 



cos — cospj, 



( ~~ 1 ) ™ \Jcos 2 p x + cos 2 p y — cos p x cos py with n = 1,2. 
The Fermi sea corresponding to the lower eigen energy 
value is given by the union of the region — 7r/2 < 
k x < 7r /2 (with all values of k y allowed) and the re- 
gion — 7r/2 < k y < tt/2 (with all values of k x allowed). 
Whereas the Fermi sea corresponding to the higher eigen 
energy value X\ is given by the intersection of the region 
— 7r/2 < k x < tt/2 and the region — n/2 < k y < tt/2. 
Since the number of electrons is equal to the number of 
sites, the total area occupied by both Fermi seas is equal 
to the area of the Brillouin zone (47r 2 ). Furthermore, the 
Fermi surface corresponds to A^ = 0. 

Using mean-field approximation we get, after averaging 
the Hamiltonian over phonon coordinates, the following 
effective Hamiltonian (see Ref. [11] for details): 



H = #1 - 2— ]T [(b\jb 2j + bl^M^ + bljhj) 
j 

+ (b\ j b lj -bt j b 2j )(b\ j b lj -bl j b 2j ) 



0J 



(2) 



where < .. > implies averaging over the phonon coor- 
dinates, At/2Mvo is the electron- JT coupling constant, 
2A(b\ j b 2j + bl j b lj ) = -w V2M^(Q 2j ), and 2A(b\ J b lj - 

bljb 2j ) = -w V2Mw (<93j>- 

We will now determine the wavevector for long range 
cooperative ordering of the Q 2 ,3 modes. For this we must 
figure out the values of q that make the susceptibilities 
X2,3(<l,q) diverge. The expressions for X2,3{o) are given 
as 



Xi(q) = -2 
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where (ci T ,ci T ) = (&!-., b\-S) ■ M, M is the diagonalizing 
matrix for the kinetic matrix T with Mi 1 = sin(^/2), 
M 2 , 2 = -sin(%/2), and M 1>2 = cos(%/2). At K 
the second term on the right hand side of both \i and 
X3 produces a divergence at Q because of the fact that 
^k+Q _ an( j t ^ at p crrm energy is zero. Fur- 

thermore the ratio of X2/X3 at K and wavevector Q is 
expected to be 3. 

3. Results and discussion 

In Fig. 1 we plot X2,3(<3) as a function of temperature 
(with hopping term t set equal to 0.1 eV) and notice that 
they diverge as T —* with X2 diverging faster than X3- 
In Fig. 2 we plot XzMQ'Q) as a function of q (with < 



q < n) at fixed temperature T = 0.01 K and hopping term 
t = 0.1 eV and find that they both peak sharply close to 
Q with X2 peaking faster. Then based on the fact that 
the phonon mode goes soft at wavevector Q we compute 
the ground state energy when only either Q 2 mode or Q3 
mode gets excited cooperatively in the system. The order 



parameters are given by {b\jb 2 j + bhbij) = c 2 cos(Q ■ Rj) 



and (bljbij - b' 2j b 2 j} = c 3 cos(Q • Rj) with -1 < c 2 ,3 < 1 
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Figure 1. Plot of X2,3(7r,7r) as a function of temperature at 
t = 0.1 eV. 
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Figure 2. Plot of X2,3(<l,<l) f° r different values of (0 < 
q < tt) at T = 0.01 K and t = 0.1 eV. 

and Rj being the position vector. Here it should be 
pointed out that the order parameter (b\jb 2 j + b\jb\j) 
corresponds to the density difference of electrons in the 
two orbitals tpx = (V^ 2 -?/ 2 — V'3z 2 -r' 2 )/v / 2 and ipy 



2 



— (ip x 2_ y 2 + ■0 3z 2_ r 2)/v / 2 (see Ref. [5]). From the sym- 
metry of the ipx an d ipY orbitals it follows that at each 
site the total charge is unity. The unit cell needed to 
compute the ground state energy consists of two adja- 
cent sites with the Brillouin zone being given by — ir < 
(k x + k y ) < 7r and — ir < (k x — k y ) < ir. We diagonalize 
a 4 x 4 matrix at each momentum and integrate the low- 
est two eigen energies over the Brillouin zone to obtain 
the ground state energy. The results of 
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Figure 3. Dependence on dimensionless polaronic energy 
(A 2 /wot) of dimensionless ground state energy per site (E/t) 
for cooperative Q2 and Qz modes. 
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Figure 4. Variation of coefficients 02,3 of ODW order param- 
eters for Q2 and Q3 distortions as a function of the dimen- 
sionless polaronic energy (A 2 /ujot). 

our calculations are shown in Figs. 3 and 4. From Fig. 
3 we see that the ground state energy corresponds to the 



Q2 mode with the difference in energy between the Q2 
only state and the Q3 only state peaking at intermediate 
values of the dimensionless polaronic energy (A 2 /ujot). 
For zero values and infinite values of the polaronic energy 
both modes yield the same energy because zero value 
implies no phononic coupling effect while infinite value 
corresponds to localized polarons. Thus for large values 
of the polaronic energy, the ground state energy is only 
slightly smaller than the polaronic energy. Furthermore, 
from Fig. 4 we also see that as the polaronic energy 
increases the values of c 2 ,3 increase and become unity 
around A 2 / (uJoi) <~ 2 implying that for the Q3 (Q2) mode 
ip x 2_ y 2 (ipx) orbital is occupied fully at one site with the 
ip 3z 2_ r 2 (ipy) orbital being fully occupied at the adjacent 
sites. 

In conclusion, we have studied orbital ordering for the 
ground state of the undoped manganite systems. We 
find that the two-dimensional orbital ordering, in the fer- 
romagnetic planes of the observed A-type antiferromag- 
netic state, is governed by the wavevector Q with Q2 JT 
mode being cooperatively excited in the system. 
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